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Abstract. We consider the vector and scalar form factors of the charm-changing current responsible for 
the semileptonic decay D — > ttIv. Using as input dispersion relations and unitarity for the moments of 
suitable heavy-light correlators evaluated with Operator Product Expansions, including 0(a 2 s ) terms in 
perturbative QCD, we constrain the shape parameters of the form factors and find exclusion regions for 
zeros on the real axis and in the complex plane. For the scalar form factor, a low energy theorem and 
phase information on the unitarity cut are also implemented to further constrain the shape parameters. 
We finally propose new analytic expressions for the Dn form factors, derive constraints on the relevant 
coefficients from unitarity and analyticity, and briefly discuss the usefulness of the new parametrizations 
for describing semileptonic data. 



1 Introduction 

The charm-changing current responsible for the semilep- 
tonic decay D — > trlv is characterized by the vector and 
scalar form factors f+(t) and fo(t), defined by 



■(p')\d^c\D°(p)) = (p'+p) M / + (t) + (p-pV-(t) ) 



fo(t) = /+(*) + 



Ml 



Ml 



(1) 



/_(*), t = q 2 = {p-v'f. (2) 



In the isospin limit the D + — > tt° form factors are obtained 
from f±(t) by multiplying with l/y/2. 

The knowledge of the shape of the Dtt form factors in 
the physical region Mf < t < (Md — M^) 2 is of interest 
for the determination of the element \V c d\ of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix entering precision tests 
of the Standard Model. Recent measurements of the branch- 
ing fractions of the semileptonic decays D — > irlv and 
D — > Klu by the CLEO collaboration [TJ[2] renewed the 
interest in the theoretical study of these processes. 

Earlier studies on the heavy-light form factors are based 
on simple pole parametrizations which implement heavy 
quark scaling laws [3J. The charm-changing form factors 
have been studied also in phenomenological models based 
on heavy quark and chiral symmetry [TJ, and more re- 
cently in the framework of Quantum Chromodynamics 
(QCD) light-cone sum rules (LCSR) [5]. Lattice studies 
have been carried out for the D — > tt, K form factors in 
the whole physical region [6j[jj[8]|9j . The D — > ttIv decay 
has also been considered recently in the emerging field of 
hard pion effective theory [TO] . 



Analyticity and unitarity are useful tools for improving 
the knowledge on the form factors. The standard disper- 
sion relations are of little use in the present case due to 
the scarce information on the Dtt form factors on the uni- 
tarity cut. However, the method of unitarity bounds [TT] 
IT2] proves to be a useful approach in cases such as this. 
The method compensates for the lack of experimental in- 
formation on the cut by an upper bound on the modu- 
lus squared of the form factor, obtained from unitarity 
and a dispersion relation for a suitable correlator of the 
same current, which can be evaluated by Operator Prod- 
uct Expansion (OPE) in the spacelike region. Employing 
standard mathematical techniques, one can then correlate 
the values of the form factor and its derivatives at differ- 
ent points inside the analyticity domain. Various versions 
of the method were applied to the pion electromagnetic 
form factor gMMIMMni, the Kir form factors QM9J 
UniinjIlllSllllllSlllS], as well as to the heavy-heavy [571 
HHUmiSnilSIlES and heavy-light form factors [SIGHIES]. 
A review of the method was presented recently in [25] . 

The D — > tt, K form factors were investigated with this 
method some time ago [33]. In this paper we revisit the 
problem, motivated in part by the progress in perturba- 
tive QCD calculations, which yield now the heavy-light 
correlators of interest to order a 2 [3jj]. Since the correla- 
tors are given in [36] only for a massless light quark, we 
restrict our study to the D — > tt form factors. 

In sect. [5] we briefly review the method of unitarity 
bounds in the heavy-light sector, extending previous stud- 
ies by exploiting also higher moments of the corresponding 
correlation functions calculated in perturbative QCD |36) . 
The simultaneous use of several independent constraints 
is expected to increase the strength of the predictions. In 
sect.lH we show also how to incorporate additional infer- 
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mation on the form factors in the analyticity domain or 
on the unitarity cut. The numerical input of the method 
is reviewed in sect. [3] 

In sect. 3] we investigate the model-independent con- 
straints imposed by analyticity and unitarity on the low 
energy behaviour of the form factors. Specifically, we con- 
sider the shape parameters entering the Taylor expansion 
around t = 0, 



Mi 



k = +, 0, 
(3) 

and derive allowed ranges for the slopes X' k and the cur- 
vatures We work with dimensionlcss parameters, the 
choice of M 2 as normalization scale being merely a con- 
vention which does not influence the results (other choices, 
for instance M|>, would simply scale the coefficients by a 
constant factor). In the same section we further apply the 
same formalism in order to isolate regions on the real t- 
axis and in the complex t-plane where zeros of the form 
factors are excluded. The knowledge of the possible ze- 
ros is of interest, for instance, for the dispersive methods 
based on phase (Omnes-type representations) and for test- 
ing specific models of the form factors. 

In sect. [5] we propose a new parametrization of the 
Dtt form factors, which generalizes the systematic expan- 
sion proposed in [3 5) for the Btt case by properly taking 
into account the position of the singularities generated by 
the first charm excited states. We derive constraints im- 
posed by analyticity and unitarity on the coefficients of 
the parametrization and demonstrate its usefulness by fit- 
ting a sample of data points generated from the CLEO 
experimental results pQ. Finally, in sect. [6] we summarize 
our conclusions. 



2 Outline of the method 

We start with the heavy- light invariant amplitudes II + (q 2 ) 
and 77o(<7 2 ) defined by the vector- vector correlation func- 
tion 

- tfgT - q»q v )n+{q 2 ) + q»q v n (q 2 ) 

= i [ d 4 xe iqa: (0\TV' J -(x)V^(0)\0) (4) 



where V u = d'y^c. It is convenient to consider the moments 
of the invariant amplitudes at q 2 = 0, 

which satisfy dispersion relations of the form 
» _ 1 [°° dt lmn k (t + ie 



t n+l 



k = +,0. (6) 



where t± = (M D ± A'U) 2 . 

From QCD it is known that the amplitude II + (q 2 ) 
satisfies a once subtracted dispersion relation, while for 



77o(q 2 ) an unsubtracted relation converges. Therefore, the 

quantities x+ and Xo are defined for n > 1 and n > 0, 
respectively. 

The connection with the form factors /+(t) and fo(t) 
defined above is provided by unitarity: including in the 
unitarity sum for the spectral functions the contribution of 
the Dtt states in the isospin limit leads to the inequalities 

Img + (t + i6)>^ [(f - f+)( *- f - )]3/a |/ + (*)| 2 , (7) 



ImFI (t + ie) > 



2 48tt t 3 

3 M- [(t-t + ){t-t^/ 2 

t 3 



\fo(t)\ 2 , (8) 



2 16tt 

which hold for t > t + . 

On the other hand, xi™' can be calculated in OPE as 
the sum of perturbative (FT) and nonperturbative (NP) 
contributions: 



In) (n)PT , (n)NP 

xl = Xk +xl ■ 



(9) 



The perturbative parts of the moments of heavy-light cor- 
relators for n < 7 were calculated up to two loops in |36j . 
Using these results we can write 



X k 



<n)PT 



(n) O) 
C fc,0 C k,l a s 



(«) 2 
C k,2 a si 



(10) 



Ob- 



where a s is the strong coupling. The coefficients c k , 
tained using eqns. (34), (35) and the Appendix of [35] are 
compiled in Table [T] for several moments that will be used 
in this work. 

The leading nonperturbative contribution of the quark 
and gluon condensates can be obtained from [541137) . and 
is written as 



n)NP 



1 



2(n+2) 



(n)NP 

xh 



(aG 



2\1 



12tt 



(n)NP 

-x+ ■ 



(11) 



(12) 



From the dispersion relations ([B]) and the unitarity con- 
ditions © and (|H]), it follows that each form factor fk(t) 
satisfies a set of integral inequalities written as 



dt Pk n) (t)\f k (t)\ 2 <xl 



(<>) 



+,0, (13) 



where the weights p k (t) are the product of l/t n+1 with 
the phase space factors entering the unitarity relations. 
We use now the fact that the form factors f+(t) and /o(*) 
are analytic functions in the complex t-plane cut along the 
real axis from t + to oo, and apply standard techniques to 
derive from ea. ([T5)) constraints on their values, in particu- 
lar on the shape parameters and on the regions in complex 
energy plane where zeros are excluded. First, the problem 
is brought to a canonical form by mapping the t-plane into 
the interior of a unit disk. This is achieved by the general 
conformal transformation 



2 (Mo) 



■\/t+ — to — sft-y. — t 



(14) 
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Table 1. Perturbative coefficients cj^j defined in (|10[) for the moments xi considered in the present work. 



i 





1 


2 


c (1) 
c (3) 


0.0024536 
0.0001690 
0.0000182 


0.0027977 
0.0002674 
0.0000342 


0.0066398 
0.0009104 
0.0001371 


c (0) 
c (1) 

c (2) 
' (1.1 


0.0126651 
0.0008179 
0.0000845 


0.0095139 
0.0013954 
0.0001860 


0.0045948 
0.0037717 
0.0006679 



Table 2. Input values in the OPE calculation of the moments. 



Quantity Values 



m*° 1.968 GeV 

m c (m c ) (1.27±g"gi) GeV 

m c (2GeV) 1.078 GeV 

(uu)(2GeV) ((-0.254 ±0.15) GeV) 3 

(aG 2 ) (7.0 ± 1.3) x nr 2 

a s (m c ) 0.39 



Table 3. OPE predictions for the vector correlators x+ ■ 



II 


(n)PT 

XV 






1 


0.0045547 


0.0002723 


0.0048270 


2 


0.0004118 


0.0000704 


0.0004821 


3 


0.0000524 


0.0000182 


0.0000706 



1 



which maps the complex t-plane cut along the real axis be written in a compact form as 
for t > t+ onto the unit disk |z| < 1 in the complex plane 
z = z(t,to), such that z(t + ,to) = 1 and 2(oo,to) = — 1. 
The real parameter to < t+ is arbitrary and denotes the 
point mapped onto the origin, z(to,to) = 0. In the new 
variable, the inequality (|13[) takes the form 



4V27T 



[t + -t(z,to))(t+-t ) 



2tt Jo 



d%£V)| 2 <x£°> 



(15) 



where the analytic functions g£ (z) are defined as 



(«), 



(16) 



Here t(z,t ) = t + - (t + - i )(l - z) 2 /(l + z) 2 is the in- 
verse of the function defined in (|14[) and u;l n ^(z) arc outer 
functions, i.e. analytic and without zeros in \z\ < 1, such 
that their modulus squared on the boundary, z = c l9 , is 
equal to pjf* (t(e l9 , to)) multiplied by the Jacobian of the 
transformation (I14j) . In our case the outer functions can 



-1/4 



3/2 



y/t+-t- + Jt + -i(z,to) 



v /t+ - t Q + y / <+ -i(z,t ) 



-(ra+4) 



(") (• \ 



t+ + y/t+-t(z,ta) 
V3(ml-ml)w { f(z) 



(17) 



y/t + -i(z,to)[y/t+-t- + y/t + -t{z,t )] 



(18) 

In the notation we omitted for simplicity the dependence 
on to of the functions g^\z), w^\z) and w n \z). 
The analytic functions g^\z) admit the expansions 



( n ) / \ 



-,(«) I „(«) 



(")^2 



(19) 



convergent in \z\ < 1. From eq. (fT5]) it follows that the 
coefficients satisfy the inequality 



E(C/) 2 < x k 

3=0 



(") 



+,0. 



(20) 
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Table 4. OPE predictions for the scalar correlators Xo 



n 


(n)PT 

Xo 


(n)JVP 

Xo 


Xo 





0.0170744 


-0.0010543 


0.0160201 


1 


0.0019357 


-0.0002723 


0.0016634 


2 


0.0002586 


-0.0000704 


0.0001883 



Since each term in the left side is positive, the largest do- 

(n) 

main allowed for the first Taylor coefficients g y k . , < j < 
J — 1 is obtained from ([20]) by setting the higher terms to 
zero. More generally, a rigorous correlation between these 
coefficients and the values = g\. (z p ) at some real 
points z p , 1 < p < P. is given by the determinantal in- 
equality [HE2III25] 



-(«) 

A n ) 

Sfc,l 



tin) 
,,2'J 



An) {ZlZ 2 ) 
Sfc 2 



An) 
?fc,2 

(ziZ 2 ) J 

1 - ZlZ 2 



1 - Zl z 2 1 



(^1^P) J (^2gp) J 

?fc,P 



?fc,P 

(zizp)* 

1 — Z\Zp 
(Z2Z P ) K 

1 - z 2 z P 



.2 ./ 



1 — ZiZp 1 — Z2-2P 



> 0. 



(21) 



where 



and 



j-i 



An) 

Sfc,p 



*(») 
Sfc,p 



(n) 



J=0 



J-1 

i=0 



p=l,2,...P. 



(22) 



(23) 



The condition (|21[) is expressed in a straightforward way 
in terms of the values fk(t P ) of the form factors at t p = 
t(z p ,to) and the derivatives at t = 0, using eqns. (|14p and 
(1161) . The generalization to complex points z p can be found 

in nng. 

In the present work we use the inequality (|21[) to ob- 
tain bounds on the slopes and curvatures of the form fac- 
tors defined in eq.©. As input we shall take the values 
/+(0) = /o(0) at t = 0, which arc known from LCSR [5]. 
An additional piece of information for the scalar form fac- 
tor is provided by a low-energy soft-pion theorem of the 
Callan-Treiman |38] type, proved in [39] , which in the Dtt 
case reads 

fo(A Dv ) = f D /U, (24) 

where Ad-k = Mp> — M% is the relevant Callan-Treiman 
point and fa and are the meson decay constants. 

As shown in |25j . from ea.([2T|) one can derive also re- 
gions where zeros of the form factors are excluded: this is 
done in a straightforward way by including in eq. (|21[) the 
input fk(zo) = and finding the values zq for which the 
inequality is violated. 

The constraints on the shape parameters can be fur- 
ther improved if some information on the form factors on 



the unitarity cut is available, in particular, if the phase 
5k(t) defined as 



h{t + ie) = \f k (t)\e l5 ^\ k = +,0, 



(25) 



is known along a low-energy interval, t + < t < ti„ . The 
implementation of this information can be done by the 
technique of generalized Lagrange multipliers and involves 
the solution of an integral equation, applied first to the Kit 
form factors in [T51[T51 . A review of the method and more 
references are given in |25j . 

Previous work on unitarity constraints for the heavy- 
light form factors 33,34, 35 , exploited dispersion relations 
for the lowest moments of the correlators, corresponding 
to n = 1 (n = 0) for the vector (scalar) form factor. 
Here we use also the higher moments calculated to two 
loops in pcrturbative QCD [36]. Specifically, we use the 
dispersion relations for the moments X+ \x+\x+ an d 

Xo°\ Xo^j Xo" 1 ■ From the inequalities (|T5l) we obtain, for 
each form factor, a family of three different constraints. 
The final allowed domain for the parameters of interest 
will be the intersection of the three individual domains. 

Up to now we did not specify the parameter io appear- 
ing in the conformal mapping (|14D . One can show that the 
bounds presented above are independent on t$: indeed, 
the bounds are obtained by solving extremum problems 
upon the class of analytic admissible functions satisfying 
an L 2 norm condition like (fP5"j) . Changing amounts to 
mapping a unit disk to another, and by this the class of 
admissible functions is not changed [40]. For deriving the 
bounds reported in sect. 2] we worked with to = 0. 



3 Choice of the input 

We have presented in Table [T] the coefficients cl". from 
[36] , entering the perturbative calculation of the moments 
■ In Table [5] we compile other quantities entering as 
input in the calculation of the moments. We exploited the 
approximate scale invariance of the product m c (uu) and 
evaluated it at a scale of 2 GeV. We took the masses and 
the strong coupling constant from the PDG tables [H] 
and used the renormalization group equations to evolve 
them to the relevant scale. The gluon condensate (aG 2 ) 
has been taken from [4"2"ll43] . The denominator of ea. (fTT|) 
involves the pole mass [37]. The PT and NP contribu- 
tions and the total OPE predictions for x^ an d Xo are 
summarised in Tables [3] and S] for the values of n specified 
above. 

As mentioned earlier, we work in the isospin limit, tak- 
ing for convenience for the D meson the mass Mp> = 
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1.869 GcV of the neutral D meson, and for tt the mass 
M w = 0.1395 GeV of the charged pion. Then t + = 4.02 GeV 2 
and f_ = 2.98 GeV 2 . 

In our analysis, we used as input the value /+(0) = 
0.67±g-J? provided by LCSR [50, and / (0) = /+(0) as 
follows from eq.@. Using f D = (206.7 ± 8.5 ± 2.5) MeV 
and U = (130.41 ± 0.03 ± 0.20) MeV gj, the low-energy 
theorem ([M]) gives f (A Dx ) = 1.58±0.07 as quoted in [SJ. 

In order to obtain the phase Sk{t) defined by (j2"5j) . a 
method useful in the case of the pion electromagnetic or 
the Kir weak form factors is to invoke Fermi- Watson the- 
orem and take the precisely known phase shifts of the cor- 
responding clastic partial waves. Since in the Dtt case the 
elastic scattering is not yet investigated (except some com- 
ments on the S- wave in [H]), we shall roughly estimate 
the phase of the form factors from the masses and widths 
of the resonances dominant at low energies. Namely, from 
the relativistic Breit-Wigner parametrization we obtain 



S(t) = arctan 



( M R r(t) \ 



(26) 



where Mr is the mass and r(t) the energy-dependent 
width 

written in terms of the angular momentum J, the width 
Fr and the cm. momentum q(t) = ' (t — t-)(< — t + )/At. 

The lowest vector D* and scalar excited states that 
couple to the Dtt system produce singularities above the 
threshold, on the second Riemann sheet. The central val- 
ues of the masses and widths of the lowest charged Dtt 
vector and scalar resonances listed in [UJ arc: Mq-> — 
2010.25 ±0.14 MeV, r D . = 96±4± 22 MeV and M D . = 
2403 ± 14 ± 35 MeV, r D » = 283 ± 24 ± 34 MeV. 




Fig. 1. Phase of the Dtt scalar form factor using a relativistic 
Breit-Wigner parametrization of the lowest resonance Dq . 



We note that the vector resonance D* is very close to 
the threshold, so that a reasonable expression of the phase 



1 Lattice simulations are in agreement with this number: 
a relativisitic computation on a fine lattice in the quenched 
approximation [B] led to the value /+ (0) = 0.74(6) (4), while the 
result /+(0) = 0.65(6) (6) was obtained in [9] using maximally 
twisted Wilson fermions with Nj — 2. 



cannot be obtained from a Breit-Wigner parametrization, 
but in the scalar case we can assume that the phase So(t) 
is reliably described by the expression (f2"6"f with J = 0, 
which is plotted in Fig. [TJ In our work we shall implement 
the phase up to the point t in = (2.6 GeV) 2 , close to the 
first inelastic Drj channel opening at 2.42 GeV. 



4 Constraints on the shape parameters and 
zeros 

The interior of the ellipses shown in Fig[2] represent the 
allowed domains in the slope-curvature plane for the vec- 

(n) 

tor form factor, obtained with three moments x+ of the 
vector correlator. We use the inequality (|2~Tj) with no other 
input on the cut or in the analyticity domain except the 
value of /+(0) given above. The best (smallest) domain 
is given by the lowest moment, but the higher moments 
contribute to slightly reducing it, since one must take the 
intersection of all the domains in order to fulfill simulta- 
neously the constraints. We also indicate the slope and 
curvature of the simple pole ansatz [3] 



/+(*) 



/+(0) 



{l-t/M D ,)(l-a Dw t/M D ,) 



(28) 



with the parameters a D7T = O^lljj;^ and M D , = 2.007 GeV 
proposed in [SJ. As seen from the figure, the point satisfies 
the unitarity constraints. 

In the scalar case we use additional information on 
the phase and the soft pion theorem given by ea. ([M|) . Fig. 
[3] illustrates the effect of these additional constraints, in 
the particular case of the bounds derived from the lowest 
moment ■ W e also show the slope and curvature from 
the pole ansatz [5] 



/o(*) 



/+(0) 



i-t/(p Dn M 2 D ,y 



(29) 



with the parameters /3dtt = 1-41 ±0.06 ±0.07 and Mo* = 
2.318 GeV suggested in [5]. The point satisfies the con- 
straint imposed by the lowest moment, taking into account 
also the phase and the low-energy theorem (|2~4"|) . 

In Fig. @] we show the allowed domains obtained by 
using different moments Xo °f the scalar correlator. The 
phase and the low-energy constraint (|24[) have been in- 
cluded in all cases. Again, the smallest ellipse is obtained 
with the lowest moment Xq ', but the simultaneous con- 
straints reduce further this domain. We obtain rather small 
allowed regions in the slope-curvature plane when the phase 
as well as the low energy theorem are simultaneously taken 
into account. In particular, the point corresponding to the 
pole ansatz (|2"9"1) falls outside the allowed regions imposed 
by the moments with n = 1 and n = 2. 

We mention that the above bounds were obtained with 
the central values of the input parameters. By varying si- 
multaneouly all the input values we can obtain more con- 
servative regions, which are slightly larger than the do- 
mains shown in the figures. The modification of the mo- 
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Fig. 2. Constraints on the slope and curvature of the vector 
form factor obtained using as input different moments of the 
correlator. The point indicates the slope and curvature of the 
pole ansatz 1)2811 . 



Ret 



■ standard bound 
standard bound + CT 
phase + CT 
pole fit 




Fig. 3. Constraints on the slope and curvature of the scalar 
form factor obtained with the moment Xo i from the standard 
bound and by including the phase and the low-energy theorem 
(1241) . The point indicates the slope and curvature of the pole 
ansatz i|2"9")) . 



3 
2 
1 

"o 



pole fit 



0.2 



0.4 0.6 

X'x 10 2 



Fig. 4. Constraints obtained by considering different moments 
of the scalar correlator, including information about the phase 
and the low-energy theorem (|24|l . The point indicates the slope 
and curvature of the pole ansatz ((29]). 



ments entering the inequalities (|15[) affects the re- 
sults in a monotonous way: larger values of lead to 
weaker bounds. In our study we reduced this source of 
uncertainty by using more precise perturbative values of 
heavy-light current correlators, calculated in [3B| at order 
a 2 . We have also studied the influence of the uncertain- 
ties in the resonance parameters for the mass and width of 
the scalar resonance |41) , and the uncertainty on the form 



Fig. 5. Domain without zeros for the vector form factor, ob- 



tained from the lowest moment x+ an d the input /+(0) 



E o 



-10 



Fig. 6. Domain without zeros for the scalar form factor, 
obtained from the lowest moment Xo an d the input /o(0) 
(smaller region) and using in addition the constraint 
(larger region). 



factor values at t = and at the Callan-Treiman point 
(fiM)) . For instance, for the lowest moment illustrated in 
Fig. [3l the union of the ellipses resulting from the vari- 
ation of j+(0), /o(Z\_D7r), Mo* and _Td» within the er- 
rors quoted in sect. |3l and by enlarging x+ by 10% yield 
the allowed ranges 0.08 • 10~ 2 < X' < 0.67 ■ 10~ 2 and 



0.51 • 10~ 5 < A[,' < 2.56 • 10~ 5 , somewhat larger than 
ranges corresponding to the central values of the input, 
0.12 ■ 10~ 2 < X' < 0.58 ■ 10~ 2 and 0.69 • lO" 5 < Aq < 
2.25- 10~ 5 . 

A similar enlargement of the ellipses given by the next 
moments in Fig. 0] makes the pole ansatz compatible with 

the allowed domain given by xi 0) - It still remains outside 
the constraint yielded by > but very close to the edge 
of the allowed domain. 

As discussed in earlier works [2U1I17]. the knowledge 
of the zeros is useful for specific fits and for testing gen- 
eral ideas about form factors. Although chiral symmetry 
for instance implies the existence of zeros in scattering 
amplitudes, there is not much information theoretically 
about them in the case of form factors. Nevertheless, phe- 
nomenological dispersive analyses often assume that zeros 
are absent in their fits to experimental data. Whereas re- 
mote zeros do not make any appreciable effect on these, 
nearby zeros if present will affect the behaviour of such 
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fits and influence the number of subtractions. Therefore it 
is important to exclude zeros in the low-energy region. 

For illustration, we consider in the following only the 
constraints imposed by the lowest moments, correspond- 
ing to n = 1 in the vector case and n = in the scalar one, 
using as input the quoted value of /+(0) from LCSR. In 
the vector case we find that zeros on the real axis are ex- 
cluded in the range (-1.0, 0.80) GeV 2 . For the scalar case 
the range with no zeros is (—2.51, 1.55) GeV 2 for the stan- 
dard bounds, while with the inclusion of the low-energy 
constraint (flM)) the range is increased to (—3.55, 3.89) GeV 2 . 
In Figs.[5]and[6]we present the regions of excluded zeros in 
the complex t-planc for the vector and scalar form factors. 
In the scalar case the zeros are excluded in a larger domain 
if the low-energy constraint (|24|) is also imposed. In princi- 
ple it would be possible to include also the phase and solve 
the integral equations for each point that is being tested. 
While an improvement of the size of the domains is ex- 
pected, the computation is a laborious one and is beyond 
the scope of the present investigation. 



5 New parametrization of the Dn form 
factors 

Until recently, the most popular parametrization of the 
heavy-light form factors was the pole ansatz [3], given in 
eqs. (f2"5f and (|2"9"|) . A more general parametrization, based 
on a systematic expansion in powers of a conformal vari- 
able, was proposed in [35] for the Bit vector form factor. 
In [5] the same type of parametrization was adopted also 
in the Dir case. Specifically, the vector Dtt form factor was 
written in [5] as 



/+(*) 



1 



1 - t/Ml, 



E^ 



(30) 



where as above z — z(t, to) is the conformal mapping (|14l> . 
Similarly, for the scalar form factor, we write 



where z — z(t, to) is the variable defined in (fT4")l . However, 
this straightforward generalization is not entirely consis- 
tent, since the pole due to the lowest D* resonance is above 
the threshold on the second Riemann sheet, while in the 
expression ([3D]) the singularity is on the real axis. Due to 
this fact, as we shall see, unitarity constraints on the free 
parameters bj cannot be derived. We mention that these 
constraints are useful not only for restricting the range of 
the independent parameters, but also for estimating the 
truncation error |35j . 

In the present work we write down improved parametriza- 
tions that take into account the proper position of the sin- 
gularities produced by the first charm excited states. We 
propose for the vector form factor, instead of eq. ((30"]) . the 
representation^ 



/+(*) 



Ml, 



J 



b+jZ 3 , (31) 

m 2 , -t + y/i- t/t+M D .r D . ^ ' J 



It is not necessary to use a P-wave phase-space in the 
Breit-Wigner (BW) denominator, since the proper behaviour 
at threshold will be imposed below for the entire function f+(t). 



fo(t) = 



Mf>, 



A/ 2 , 



t + y/i-t/t+M D .r DS 



3=0 



The coefficients bk,j are the parameters to be used in fits 
of data. Unitarity and analyticity imply that they are not 
completely free, but satisfy a constraint |35j . In order to 
derive it, we insert the representations (|5Tj) and (|32p into 

the definition ([TB"]) of the functions g^\z), extract the 
corresponding Taylor coefficients and apply the inequality 
(f2T)]) . For simplicity, we consider only the constraints ob- 
tained with the lowest moments of the correlators. Then 
the above steps lead to the inequalities 



k = +,0, 



(33) 



where B^f* arc calculated as 



B 



a-) 



E Vk,mVk,m+\i-j\ , k = +, 0. (34) 



Xq ' , and the 



Here we used the notations x+ = X+\ Xo 
numbers rjk, m are the Taylor coefficients appearing in the 
expansions 



Ml,w^(z) 



M 



D> 



t(z,t ) + y/l-i{z,t Q )/t+M D .r L 



E 



(35) 



M 



' 2 , - t(z, to) + y/1 -t{z,to)/t+M D *r E 



2^ ^ 

m=0 



(36) 

where the outer functions w^\z) are defined in eqs. ([T7| 
and (fT5|) . and t(z,t n ) is the inverse of (fT4"l) . 

As remarked in sect. [2 the bounds derived in the previ- 
ous section are independent of the choice of the parameter 
to (in the calculations we have set t to 0). By contrast, 
the parametrizations given above are based on truncated 
expansions and depend on to- For to = the scmilcptonic 
region (0,£_) is mapped onto the range (0,0.35) in the 
z-plane. This may lead to quite large truncation errors at 
the right end of the physical region. As suggested first in 
30 , 33] , it is more convenient to choose to such as to map 
the physical range symmetrically around the origin of the 
z-planc, i.e. z(0,to) = —z(t-,to). This gives the optimal 
value 



opt 



1.97 GcV^ 



(37) 



when the physical range is mapped onto the segment (- 
0.17, 0.17) in the z-plane. In the following we shall inves- 
tigate both choices t = and t n = t opt . 

Since the outer functions and the BW representations 
are analytic in the unit disk and have no singularities on 
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the boundary \z\ = 1 (the denominators vanish outside 
the unit disk), the expansions in the right hand sides of 
eons. ([35"]) and ([55)1 are rapidly convergent. Therefore the 

(k) 

quantities B\- can be calculated with precision by trun- 
cating the series in (|34"]) at a finite order0 

The numbers Bffl and B^' up to J = 5 for tn = are 

= 4.45 • 1CT 2 , = 2.53 • 1CT 2 , 



B, 



(+) 



02 

( + ) 
04 



-4.53 • 1CT 4 , B^ ] 



-3.60-10" 



-4.60 -10" 5 , B^P = —5.81 • 10~ 4 , (38) 



5$ =2.83, fig? =1.18, B$ = -1.68 



B 



(o) 

03 



-2.29, B$ = -0.32, B^> = 1.60. (39) 



? (0) 



while for to = t opt the values read 



B$ = 4.45 • 10" 



B { + } = 1.69 • 10" 



B, 



( + ) 

02 



-8.51 ■ 10" 



B, 



(+) 

03 



-8.87- 10~ 4 , 



B J> = 5.75 • 10 



B { + } = -3.10- 10~ 3 , 



(40) 



,(0) 
00 
(0) 
03 



1.91, B, 



(o) _ 

L — 

>(0) . 



0.25, B 



(o) = 

(0) 



-1.61 



-0.53, B^ - 1.09, B 05 



0.51. (41) 



The remaining coefficients are obtained from the rela- 
tions B$ + js = Bq) and the symmetry property B^ = 



B^ ' which follow from ([34"]) . In the calculation we used for 
the masses and widths of the resonances the central values 
from [41], quoted above. The smaller values of B^ are 
due partly to the value of the QCD vector correlator and 
the form of the unitarity inequality ([7|l. and partly to the 
parameters of the D* -resonance. 

Additional information, like the values of the form fac- 
tors at t = and the low energy theorem for the scalar 
form factor can be implemented easily. For the vector form 
factor it is useful to implement also the P-wave type be- 
haviour at the unitarity threshold t+ . As discussed in [35] , 
this condition is equivalent to 



d/+(t(Mo)) 
dz 



0. 



and amounts to the simple algebraic relation 



(42) 



(43) 



7=0 



3=0 



which must be fulfilled by the parameters fo+.j, simultane- 
ously with the condition (|3U)) . These constraints are useful 
in the fits of the data, reducing the number of independent 
parameters. 



3 A pole situated on the real axis, as in (|30|) . is mapped by 
(|14|l onto a point Zo on the boundary of the unit disk, |zo| = 1- 
The quantities rjk, m involve the powers 1/z™, and it is easy to 
see that the series in the right side of (|34[) is not convergent in 
that case. 




Fig. 7. Fit of a sample of data on the vector form factor /+ (t) 
in the semileptonic region with the parametrization (|31|) . 



In order to illustrate the properties of the new parametriza- 
tions, we generated from the CLEO data [1] a sample of 
values and errors for the vector form factor f+(t). We 
obtained the values from Table XIII of [T], using for con- 
venience \V c d\ = 0.236 as in pp. Assuming that the quoted 
values correspond to the center of the bins defined in Ta- 
ble IX of PP, and adding /+(0) from LCSR as quoted in 
sect. [31 we obtained the 10 data points given in Table [SJ 
We then made a fit of these data using the expressions 
([3T|l . with the threshold condition (|43[) . using both to = 



and to = t opt . The results of the two choices are actu- 
ally very similar: with 2 independent parameters 6 + .o an d 
the third one 6 +i 2 being determined from ea.([43)). we 
obtained for to = a minimal x 2 = 9-71 with the optimal 
parameters: 

b+,o = 0.665, b+s = 0.725, b +a = -1.066, (44) 



while for to = t opt from (j37|) we obtained x = 9.67 and 
the optimal coefficients 



0.757, 6+i = 0.397, 6 4 



-0.852. 



(45) 



We have checked that the coefficients (pH)) and sat- 
isfy the constraint ([33]) with B^ } from ([38]) and ([40]) . 
respectively. The above parametrizations lead actually to 
almost identical form factors in the semileptonic region, 
as shown in Fig. [7] However, the choice to = t opt allows a 
more accurate determination of the systematic error: as- 
suming, as in [35] . that a reasonable definition of the trun- 
cation error is given by |6+, 3 | ma2: |z(i, t )| 3 , where \b +i3 \ max 
is the maximum value of the next coefficient allowed by 
the condition (|33[) at the optimal point, we find for to = 0, 
when Ifr+^lmaz = : 5.26 and z(t-,0) = 0.35, that this error 
reaches a value of about 7%, while for to = t opt , when 
|6+, 3 |maa; = 4.28 and z(t-,t opt ) = 0.17, the truncation er- 
ror can be reduced to less than 1% in the whole semilep- 
tonic region. 

The nontrivial role played by the condition Q33p is seen 
if one makes a fit with one more independent parameter in 
([3T]) . Then, for to = topt, the best fit without constraints 
gives x 2 = 6.86. However, the unitarity condition ([33]) is 
violated by the optimal parameters, and should be im- 
posed explicitly for a suitable fit. 
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Table 5. Sample of data on f+(t) used for a fit with the parametrization (|31[) . 

n 1234 56789 10 

t„ 0.15 0.45 0.745 1.04 1.34 1.59 1.89 2.23 2.68 

/+(*») 0.67 0.68 0.74 0.760 0.94 0.98 1.06 1.57 1.69 1.94 

Sf+(t n ) 0.10 0.03 0.04 0.05 0.05 0.06 0.08 0.11 0.23 0.29 



We finally remark that although the parametrization 
(|31[) for to = topt is defined such as to optimize the de- 
scription of the scmilcptonic region, it has a reasonable 
behaviour also outside this region and even on the unitar- 
ity cut. In particular, the polynomial in the numerator of 
([51]) becomes complex for t > t+, when \z\ = 1, but one 
can check that the phase of the product in (]3"TT) increases 
smoothly above the unitarity threshold and the resonance 
shape is not distorted, especially when the condition (|43l) 
is imposed. 

The above exercise demonstrates the usefulness of the 
new parametrization (|31|) for data analysis. Of course, in 
a complete study the new expressions should be used di- 
rectly in the fit of the partial branching fractions over bins 
in t = q 2 , given in [I] . The values provided by lattice calcu- 
lations in the physical region [BHIHj can be also included 
in a combined fit, allowing a simultaneous determination 
of the vector form factor and of | V c d | , as was done for the 
similar Bit case in [33] , 

6 Conclusions 

In this work, we have explored the implications of ana- 
lyticity and unitarity for the Dtt form factors. The work 
was motivated by the recent experimental measurements 
of the scmilcptonic D — > trlv decay [HE] of interest for the 
extraction of the element \V c d\ of the CKM matrix. We 
have applied the formalism of unitarity bounds, which use 
as input the dispersion relations satisfied by the moments 
of suitable heavy- light correlators at q 2 = 0. The avail- 
able calculations of the moments in perturbative QCD up 
to 0(a 2 ) terms [36] allow more precise predictions in this 
framework. We have derived a family of constraints to be 
satisfied by the shape parameters of both the scalar and 
vector form factors at t = and explored the region on 
the real axis and in the complex plane where the form 
factors cannot have zeros. The theoretical input at the 
Callan-Trciman point and the phase information for the 
scalar form factor play an important role in constrain- 
ing the corresponding shape parameters. The exclusion 
regions for the zeros that we have isolated basically cover 
a significant part of the entire low energy region. 

We have also proposed an improved parametrization of 
the Dtt form factors in the semileptonic region, by prop- 
erly implementing the singularities related to the lowest 
charmed resonances. The new parametrizations, given in 
eqns. (]3"T]) and (]3"2j) for the vector and scalar case, respec- 
tively, are based on truncated expansions in powers of the 
conformal mapping z = S(t, to) defined in (Tl4]) . with coef- 
ficients satisfying the quadratic condition (]3"3"]) . Additional 
low-energy constraints or the threshold behaviour for the 



vector form factor can be easily implemented, as shown 
above. Using a sample of values produced from the CLEO 
experimental data [T], we demonstrated the usefulness of 
the new parametrization for the description of the semilep- 
tonic data and a reasonable estimate of the systematic er- 
ror. A more complex analysis, allowing an accurate simul- 
taneous prediction of the form factor and of \V c d\, implies 
combined fits of the partial branching fractions over bins 
and of the LCSR and lattice predictions. This analysis will 
be reported in a future work. 
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